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Abstract 

We present ab-initio calculations of the magnetic moments and magnetic anisotropy energies of small FeCo clusters of 
varying composition on top of a Cu(100) substrate. Three different cluster layouts have been considered, namely 2x2, 
3x3 and cross-like pentamer clusters. The ratio of Co atoms with respect to the total number in a chosen cluster ("con- 
centration") was varied and all possible arrangements of the atomic species were taken into account. Calculations have 
been performed fully relativistic using the embedded-cluster technique in conjunction with the screened Korringa-Kohn- 
Rostoker method and the magnetocrystalline anisotropy energy (MAE) has been evaluated by means of the magnetic 
force theorem. A central result of the investigations is, that the size of the magnetic moments of the individual Fe and 
Co atoms and their contributions to the anisotropy energy depend on the position they occupy in a particular cluster 
and on the type and the number of nearest-neighbors. The MAE for the 2x2 and 3x3 clusters varies with respect to the 
"concentration" of Co atoms in the same manner as the corresponding monolayer case, whereas the pentamer clusters 
show a slightly different behavior. Furthermore, for the clusters with an easy axis along a direction in the surface plane, 
the MAE shows a significant angular dependence. 

PACS numbers: 73.20.At, 72.10.Fk, 73.22.-f, 75.30.Hx, 73.20.Hb 
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I. INTRODUCTION 



Surface-supported nanoparticles of magnetic atoms arc generating a lot of interest nowadays, because of 
potential applications in non-volatile magnetic storage media. In the last few years, magnetic nanostructures 
were investigated experimentally in terms of various different methods like STM, XMCD and MOKEi^££. In 
using these techniques together with phenomenological models^ and sum-rules^ it was possible to infer the high 
anisotropies and orbital moments of single magnetic adatoms on a non-magnetic substrate. Eventually it is the 
ambition of experimental methods to produce and manipulate nanostructures on an atom-by-atom level which 
is currently attempted by means of for example magnetic tunnel tips. 9 Furthermore, by combining magnetic 
and non-magnetic materials like CoPt 1 , or two different magnetic species like FeCoiS to form nanoclusters, 
interesting structural arrangements and enhanced, tunable magnetic properties are obtained. 

A key challenge for experimental methods is, in view of the fabrication of so-called bit-patterned media ; 11 ' 12 
to produce grains of controlled size and position, of known composition, and with sharply defined magnetic 
propertiesj 13 ' This would path the way to increase achievable areal densities^ of magnetic recording media 
by several orders of magnitude . 13 i 16 i 17 i 18 i 19 

One of the main issues is the magnetic anisotropy energy which determines the orientation of the magnetization 
of a cluster with resect to the surface. Large MAE barriers can stabilize the magnetization direction in the 
cluster and a stable magnetic bit can be created* 2 ^! In contrast to bulk solids, clusters deposited on surfaces 
offer additional degrees of freedom to change the MAE. It can be influenced by the shape, size and composition 
of the cluster and by the substrate. In this paper we will concentrate on the investigation of the effects of 
composition and of the details of the atomic arrangements within magnetic clusters on the anisotropy by means 
of ab-initio calculations. 

The information learned from these computational studies can be used later to construct nanostructures with 
optimized properties. Moreover, these studies can be compared to calculations made for monolayers using the 
Coherent Potential Approximation (CPA) and relate the impact of the local environment to the mean field 
solution. 



II. THEORETICAL APPROACH AND COMPUTATIONAL DETAILS 



Self consistent, relativistic calculations for FeCo clusters on a Cu(100) surface have been performed using the 
embedded-cluster technique^ 2 , within multiple scattering theory (MST) which enables the treatment of a finite 
cluster of impurities embedded into a two-dimensional translational invariant semi-infinite host. Within MST 
the electronic structure of a cluster of embedded atoms is described by the so-called scattering path operator 
(SPO) matrix given by the following Dyson equation^ 2 - 

r c (e)=r h (e)[l-(t- 1 (e)~t c 1 (e))r h (e)}- 1 , (1) 

where Tc(e) comprises the SPO for all sites of a given finite cluster C embedded in a host system, th(e) 
and Th(e) denote the single-site scattering matrix and the SPO of the unperturbed host sites in cluster C, 
respectively, while tc(e) stands for the single-site scattering matrix of the impurity atoms. Once Tc(e) is known 
all corresponding local quantities, i.e., charge and magnetization densities, spin and orbital moments, as well as 
the total energy can be calculated. In all cases the atomic sphere approximation (ASA) was applied. 

Self consistency is achieved by varying the effective potentials and exchange fields using the local density 
functionals of Ceperley- Alder (in the parameterization due to Perdew and Zunger)^ 3 . and solving the Poisson 
equation as described in Ref. l22|. In all self consistent calculations for the embedded clusters the orientation of 
the magnetization was chosen to point uniformly along the surface normal [z axis). For the calculation of the 
i-matrices and for the multipole expansion of the charge densities (needed to evaluate the Madelung potentials), 
a cutoff for the angular momentum expansion of Z max = 2 was used. 

The host and the cluster sites refer to the positions of an ideal fee lattice with the experimental lattice 
constant of Cu (a = 3.6147 A). Structural relaxations of both the cluster-substrate distance as well as of the 
bond length between cluster atoms, which may also affect the magnetic properties ) 24 ' 25 ' 26 ' 27 have been neglected. 
As the atomic radii of the elemental Fe, Co, and Cu are similar, relaxations can be expected to be relatively 
small. For the semi-infinite Cu(100) host a self consistent fully relativistic calculation was performed in terms 
of the screened Korringa-Kohn-Rostoker method^ using 66 k\\ -points for the irreducible Brillouin zone (IBZ) 
integration and 16 energy points for the energy integrations along a semicircular contour in the complex energy 
plane by means of a Gaussian quadrature. 

The magnetic anisotropy energy (MAE), is defined as the difference in the total energies of the clusters with 
the magnetization along two different directions. While in the bulk or thin films of ferromagnetic materials 
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the preferred magnetization direction will be along one of the principle crystal axes, this will not be the case 
in composite nanoclusters. Depending on the shape and the composition of the clusters the easy (or hard) 
axis can be found along a (arbitrary) direction a = (sin 9 cos ip, sin 9 shop, cos 9). In terms of the magnetic 
force theorem^ the MAE is given by the energy difference between the band energies corresponding to two 
orientations, a and /3, of the magnetization 



which were evaluated using the self consistent potentials taken from the calculation where the orientation of 
the magnetization was chosen to be perpendicular to the surface (z axis). Here the band energy is actually the 
grand canonical potential, which at T — OK is given by: 



where eb is the bottom of the valence band and n(e, a) denotes the density of states for the magnetization 
pointing uniformly along direction a. Use of the force theorem in the present context seems justified as the atomic 
species considered are 3d transition metals for which the approximations due to the first order perturbation 
theory are legitimate ! 24 ' 30 ! 31 

For the present calculations only the nearest-neighbor shell of host and vacuum sites surrounding the Fe and 
Co atoms forming the clusters was taken into account. Our previous convergence studies^ 2 - of the magnetic 
properties of pure Fe clusters as a function of the number of shells demonstrated, that for a Ag substrate 
considering only the first shell suffices. For a Cu host, being a noble metal with similar properties as Ag, a first 
shell of neighbors can therefore safely be assumed to be adequate. In the following, the term 'cluster' will be 
used (for simplicity) to denote the set of Co and/or Fe atoms only. 

The dimensions of the clusters were varied in the following way: we chose a 2x2, a cross-like pentamer and a 
3x3 cluster as illustrated, e.g., in Ref. [H. The atoms in the clusters occupy nearest neighbor positions of the 
ideal fee (100) geometry. Taking the first shell, the studied systems actually consist of 30 ASA spheres in the 
2x2 case, 40 for the pentamer cluster and 53 in the case of the 3x3 cluster. These systems refer to the first 3 
layers of a parent fee lattice. The first layer is formed by atoms of the substrate, the second layer contains Co 
and/or Fe atoms and surrounding vacuum spheres and the third layer is formed only by vacuum spheres. We 
varied the 'concentration' of the Co atoms, namely the ratio of Co atoms with respect to the total number of Co 
and Fe atoms in the clusters and considered, in addition, for each concentration also all possible configurations, 
i.e. arrangements of the atoms in the cluster. 

In the following sections, the most important magnetic properties of the clusters (detailed for each constituent 
atom) are discussed, namely the magnetic spin and orbital moments and the anisotropy energy. All quantities 
have been studied with respect to the clusters' size, to the constituent atoms concentrations and with respect 
to the different positions of the atoms within the clusters. The results are furthermore compared to those of a 
monolayer of either Co or Fe, in particular cases also to a monolayer of Co^Fei^a, on top of Cu(100). 

III. RESULTS AND DISCUSSION 
A. Magnetic moments 

One of the central results of all previous studies of (pure) nanoclusters has been that the spin-only magnetic 
moments show a negligible dependence on the direction of the magnetization, while the orbital moments depend 
significantly on it . 22 ' 31 ' 34 ' 35 However, it will be shown that in clusters of composite materials, there appears 
to be a significant correlation between the size of the spin and the orbital magnetic moments and the specific 
positions of the atoms. As far as the spin moments are concerned, it will be pointed out that for some clusters 
this dependency is also connected with the atomic coordination, and hence it is due to the position occupied by 
the atom. In order to present the large amount of data in a compact way, below, the results for the magnetic 
moments are summarized separately for each type of cluster that has been studied. 

1. 2x2 cluster 

We start our discussion with the magnetic moments for the smallest cluster studied, namely for the tetramer. 
In Fig. [TJ the variation of the spin (S) and orbital moments (L) of a selected Co and an Fe atom is plotted for 
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FIG. 1: Spin and orbital magnetic moments of selected Co (left panel) and Fe (right panel) atoms as a function of 
all cluster configurations possible (Co: iA, Fe: iB, i — 1, . . . , 5), see the scheme to the right. Dashed horizontal lines 
refer to the corresponding Fe/Cu(100) or Co/Cu(100) monolayer values for a magnetization along the coordinate z-axis. 
Dash-dotted horizontal lines to those with the magnetization along the x(j/)-axis. The investigated atom is indicated by 
boldface letters. 



different cluster configurations (full lines), together with the corresponding Fe or Co monolayer values (dashed 
and dash-dotted lines). The spin magnetic moment of the chosen Co atom, even in a small 4-atom cluster, 
has almost the same value as in the Co monolayer case, whereas the Fe magnetic moments in the cluster are 
enhanced with respect to the Fe monolayer due to the reduced coordination number. In the 2x2 cluster, all 
the positions are geometrically equivalent, however, the orbital magnetic moments are still sensitive to the 
type of the atoms occupying the different positions in the cluster, or, shortly to the local environment (see 
configurations 4A and 5A, 4B and 5B). As can be noticed from the left panel of Fig. [1] the orbital moment 
values of the investigated atoms are more sensitive to the environment than those of the spin moment. The 
value of the orbital moment on an Fe or Co atom can vary by more than 30% depending on its surrounding 
atoms. While the spin moments have a very small variation with respect to the changes in the direction of the 
magnetization (from the in-plane to out-of-plane), for the orbital moments the difference can be significant. For 
the Co atoms the values of L x are higher than those for L z so it has to be expected that the Co atoms have a 
preferred in-plane direction for the magnetization. The value with respect to L z of the Co atom has a minimum 
when it has one Fe nearest-neighbor (NN) and increases with the number of Fe atoms in the cluster. Contrary, 
the same quantity of an Fe atom shows only minor oscillations when the number of the Co NN changes. 



4 



3.2 

3.1 

I 3.0 

^2.9 

2.8 

2.7 

2.6 
0.25 

5 0.20 k 



o 

i M 0.15 



0.10- 



t — i — i — i — i — rp — i — i — i — i — i — i — i — i — i — i — r 



O 

O Fe O 
O 



Fe 
O O O 
O 




J I I I I I I I I I I I I I I I I L 




ML value z 

ML value x,v 





n n< I I I I I I LU I I I I I I I I I I L 

la 2a 3a 4a 5a 6a lc 2c 3c 4c 5c 6c 7c 8c 9c 10c 1 lc 12c 



configuration 



configurations 


Fe 


Fe 




Fe 




Fe Fe Fe 


Fe Fe Fe 


Fe 


Co 


Fe 


Fe 


Fe 




Co 




la 


lc 




7c 




Fe 


Fe 




Fe 




Fe Fe Fe 


Fe Fe Fe 


Fe 


Co 


Fe 


Co 


Co 




Co 




2a 


2c 




8c 




Co 


re 




i7.. 
le 




CoFe Fe 


Fe Fe Fe 


Co 


Fe 


Co 


Fe 


Co 




Co 




3a 


3c 




9c 




Fe 


Fe 




17.. 

le 




Co Fe Co 


Pp (~*rs Pp 


Fe 


Co 


Co 


Fe 


Fe 




Co 




4a 


4c 




10c 




Co 


17.. 

le 




1.',, 

r e 




Co Fe Co 


Co Co Fe 


Co 


Co 


Co 


Fe 


Fe 




Fe 




5a 


5c 




11c 






y 










Co 


Fe 




Fe 






Co Fe Co 


Co Fe Co 


Co Co 


Co 




Co 


Fe 




Co 




6a x 


6c 




12c 



FIG. 2: Spin and orbital magnetic moments of a selected Fe atom as a function of its position in the cluster and of the 
number of Co neighbors. For the actual configuration see the panel to the right. Dashed horizontal lines refer to the 
corresponding Fe/Cu(100) monolayer values with the magnetization along the coordinate z-axis. Dash-dotted horizontal 
lines to those with the magnetization along the x(y)-axis. 



2. Pentamer cluster 



The distinctive feature of a pentamer cluster - in comparison with a quadratic 2x2 cluster - is the existence 
of two geometrically nonequivalent positions. While the atom at the central position in the cluster has four 
coordinated atoms, an atom at a corner position has only one. Fig. [5] illustrates the variation of S and L of 
particular Fe atoms with respect to changes of neighboring atoms. In the left panel the magnetic moments of 
the central Fe atom are displayed and in the right panel those of an Fe atom in a corner position are shown. 

Although the value of the spin and the orbital moments of the central Fe atoms increase monotonously with 
the number of Co NNs, they are still in a range of approximately ±0.1/j,b with respect to the monolayer value of 
2.85 /is of Fc/Cu(100). The spin moment shows linear dependence on the number of the Co NNs. In contrast, 
the values of the corner atoms, while being much higher than in the monolayer case, exhibit a less pronounced 
dependence on the environment. Taking a closer look at the orbital moments, we see that they are significantly 
larger than in the monolayer (0.078/is) only for corner atoms for which the value ranges between 0.2 and 
0.23/is, while for the central atoms they vary smoothly and are close to the monolayer value. The difference 
in the orbital moments values, namely between L z and L x is about 0.06/is and suggests that the corner Fe 
atoms prefer an out-of-plane direction for the magnetization. Finally, the spin-moments of an Fe atom at a 
corner position is increased dramatically by approximately 0.4/is as compared to the value in an Fe monolayer 
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FIG. 3: Spin and orbital magnetic moments of a selected Co atom as a function of the position occupied in the cluster 
and of the number of Fe neighbors. For the actual configuration, see the panel to the right. Dashed horizontal lines refer 
to the corresponding Co/Cu(100) monolayer values. Dash-dotted horizontal lines to those with the magnetization along 
the x(y)-axis. 



on Cu(100) {2.77 hb)- 

Fig. [3] shows the same quantities and additionally the orbital moment for the magnetization pointing long 
the y-axis, L y , now for a Co atom in the central and at a corner position. Contrary to the case of Fe, the spin- 
moment of a central Co atom is smaller than the corresponding monolayer value (1.84/ig) in every configuration. 
Similarly to the Fe case the dependence of S of the central Co atom on the number of Fe NNs is linear but with 
opposite sign. The decreased spin-moment attains its smallest value when the Co atom is surrounded by Fe 
atoms only (configuration 6a in Fig. [3j . Then again, the magnetic moment of a Co atom positioned at a corner 
is larger than in the monolayer. This is not primarily due to the enhancement of the spin- moment but to that of 
the orbital moment. In the lower panel of Fig. [3] the orbital moments are plotted for magnetizations along the 
z-direction and along the in-plane x and y directions. It can be seen, that a central Co atom shows only a small 
anisotropy of L and a slight dependence on the cluster configuration. At a corner position, a Co atom exhibits 
a much more pronounced anisotropy and for configurations lc-3c L x is largest, assuming a given orientation 
of the clusters with respect to in-plane x-y coordinate system as illustrated in the legend of Fig. [3] For all 
other configurations the orbital moment along the y-direction becomes largest and is significantly enhanced in 
comparison to the monolayer case. The reason for the strong increase of the orbital moment of the Fe and Co 
atoms at corner positions has been discussed^ for pure Fe, Co, and Ni clusters and can be traced back to the 
anisotropic environment and reduced coordination number. In the present study this effect is combined with 
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FIG. 4: Spin and orbital magnetic moments of a selected Fe atom in a 3 x 3 cluster as a function of its position and the 
number of Co NN for selected cluster configurations (see the panel to the right). Dashed horizontal lines refer to the 
corresponding Fe/Cu(100) monolayer values with the magnetization along the coordinate z-axis. Dash-dotted horizontal 
lines to those with the magnetization along the x(y)-itxis. 



the changes of the magnetic properties due to the different chemical order of the two constituents forming the 
cluster. 



3. 3x3 cluster 



The magnetic moments for the three geometrically inequivalent positions (center, edge, corner) of Fe atoms 
are displayed in Fig. H] for selected cluster configurations. Starting in each case from a pure Fe cluster, the NN 
positions are successively substituted with Co atoms until all NN sites are occupied exclusively by Co atoms. 
The only exception is configuration 2c where a 2nd NN has been considered. 

As a first observation it can be seen, that the spin-moments in the pure clusters depend on the position of the 
Fe atom. The values differ by approximately 0.1/ib when comparing them at center and center of edge positions. 
At a corner position the spin-moment is further enhanced by compared to the edge position. In addition 

there appears to be an almost linear dependence on the number of coordinated Co atoms at NN sites. While it 
is interesting to observe that the spin-moment of a central Fe atom in a pure Fe cluster is practically identical 
to the monolayer value, the largest moment of 3.06/ib is found when a Fe corner atom is coordinated by two 
Co atoms. 

The orbital moments of a central Fe atom are almost identical to their values in a monolayer, with L x being 
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FIG. 5: Spin and orbital magnetic moments of a selected Co atom in a 3 x 3 cluster as a function of its position and 
the number of Fe NN for selected cluster configurations (see the panel to the right) . Dashed horizontal lines refer to the 
corresponding Co/Cu(100) monolayer values with the magnetization along the coordinate z-axis. Dash-dotted horizontal 
lines to those with the magnetization along the x(y)-axis. 



slightly larger than L z , exhibiting negligible dependence on the local environment. Similarly, at an edge the 
orbital moments do not differ significantly from the monolayer values, and are practically independent of the 
number and position of Co NN. However, L z is now augmented in comparison to L x . Most notably the Fe atom 
located at a corner shows a much more significant enhancement of L z and is further increased as more Co NN 
are added. This can be directly related to the anisotropic surrounding at that position and reduced coordination 
number, as was already seen in the case of smaller clusters. L y of Fe atoms does not differ significantly from L x 
and therefore has not been included in Figs. [HO and [5J 

Independent of the position assumed by the Co atom in the cluster, the spin moments decrease almost 
monotonously with increasing number of Fe NN (see Fig. [5]). Except for configurations la and 2a, all values are 
smaller than the monolayer ones. 

Again the orbital moments show a dependence on the position that the Co atom occupies in the cluster. 
However, L x , L y , and L z depend in different ways on the addition of Fe atoms on the NN sites. While at 
an edge position L z is practically unchanged, L x shows an oscillatory behavior and L y experiences a strong 
enhancement as illustrated in the lower panel of Fig. [5] Then again, a Co atom at a corner has a larger L x 
almost identical to L y , and the difference to L z is approximately constant. 

In order to summarize the data discussed above, in Fig. \E\ the variation of the spin and orbital magnetic 
moments averaged over all the configurations for a given Co concentration has been plotted. One can see that 
the difference in the averaged values of the 3x3 and 2x2 clusters is very small for both S z and L z . The atoms 



S 



t 1 1 1 1 1 r 




I i i i i i i i i i i i I 

0.2 0.4 0.6 0.8 1 
Co concentration 

FIG. 6: Spin and orbital moments averaged over all the configurations for a given Co concentration. (•: 2x2 cluster, 0: 
pentamer cluster, A: 3x3 cluster) The dashed line refers to the case of a Co^Fei-a; monolayer on Cu(100). 

within a pentamer cluster exhibit a larger orbital moment with respect to the value of a monolayer of Co^Fei-a; 
on Cu(100), and also compared to the values in the 2x2 and 3x3 cases. The variation of the spin moments 
with respect to the Co concentration is linear, just as in the monolayer case, whereas the orbital moment has 
a more complex behavior. This can be related to the reduced symmetry of the local environment and to the 
higher ratio of corner atoms in these clusters. 

B. Magnetic anisotropy energy 

The magnetic anisotropy arises from that part of the total energy which is anisotropic with respect to two 
different magnetization directions^ 2 One contribution comes from the shape anisotropy (also referred to as 
demagnetization or dipole-dipole energy) which is, however, quite small (in the /LteV range) if the clusters are 
small and of square shape. Due to the strong anisotropy of the orbital moments, and the dependence of both, 
the spin and the orbital moments on the local environment, this energy is nevertheless enhanced in comparison 
to pure clusters. However, it is negligible compared to the most important contribution to the magnetic 
anisotropy in this case, namely the magnetocrystallinc anisotropy energy which is caused by spin-orbit coupling 
of the electrons in a magnetic material. This energy is rather sensitive to the local environment and due to 
the arrangement of atoms, a magnetization along certain orientations is energetically preferred. Subsequently 
it will be shown, that the magnetocrystalline anisotropy is closely related to the structure and symmetry of the 
investigated clusters. 

The MAE of the three cluster types has been investigated with respect to perpendicular orientations of the 
magnetization along the coordinate axes, as well as according to the dependence on the azimuth angle ip. To 
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FIG. 7: Magnetic anisotropy energy as a function of the Co concentration: distribution of the anisotropy energies for 
3x3 clusters as a function of the clusters configuration (diamonds); the anisotropy energy (AE±) — {(AE XZ + AE yz ) /2), 
where {) refers to the average over all configurations (solid line); anisotropy energy for Co a; Fei_ a ;/Cu(100) monolayer 
(dashed line). The index _L indicates that the MAE is calculated as the energy difference between the (surface) normal 
z and the in-plane x or y direction. 



proceed systematically, our discussion will start with the magnetization along x,y, and z, and will then focus on 
the angular dependence of the MAE. 

We compare the results for the 3x3 clusters with those for one monolayer, in Fig. [TJ where the configuration 
averaged MAE with respect to in-plane (along the i-axis) and out-of-plane (along the z-axis) magnetization, 
(AE)± = (£* — E~), is displayed as a function of the Co concentration. The size of the anisotropy energy 
is strongly influenced by specific cluster configurations. In some cases the type of the configuration may even 
be responsible for a change in the preferred orientation of the easy axis as is the case of x = 0.44% Co and 
x = 0.55% Co. Starting from a pure Fe cluster, which has a large anisotropy energy (+0.6 meV/atom), favouring 
out-of-plane orientation, the increase in the Co concentration is responsible for a decreasing anisotropy energy 
until the easy axis prefers an in-plane orientation. Then a further increase in the Co concentration leads to an 
augmentation of the in-plane anisotropy. The general trend, however, is very similar to the one obtained in a 
CPA calculation for a monolayer of Co K Fei_ x /Cu(100). In comparison to the MAE calculated for a monolayer 
of Co a; Fe(i_ :c )/Cu(100) (dashed line) the values for the 3x3 clusters are about three times larger. 

It has been suggested by Brunch and Gambardella 4, ~ that the direction of the easy axis is strongly related to 
the magnitude of the orbital moments. We investigated the orientation of the MAE with respect to the variation 
of the orbital moments for the 3x3 clusters, by choosing a configuration with high planar symmetry and applying 
it to two different concentrations (Fe4Co5, FesCo/i). In order to investigate the correspondence between the 
orbital moment anisotropy and the MAE, in Fig. [5J the values of the MAE and of the orbital moment per atom 
are plotted with respect to the positions occupied in the cluster. For the two selected clusters, three different 
magnetization directions (M||z, M||ir, M||j/ ) were considered. The rule is strictly followed in the case of FeziCos 
and stands for almost all atoms in Fe5Co4, except for the central Fe (i.e., Fej) which actually has a very small 
orbital moment anisotropy (Fig.[5J upper panel). The equivalent positions within a 3x3 cluster (labeled by the 
same number) are schematically represented in the upper part of Fig. [HJ For sake of simplicity, in this figure 
only the values for noncquivalcnt positions are presented. The nonequivalence between positions 2 and 3 is 
only due to the planar symmetry, which can also be seen from the interchange of the values for L x and L y . 
While the Fe atoms give an out-of-plane contribution to the MAE, the Co atoms have a tendency to an in-plane 
magnetization. The values of the atomic contributions (from both Fe and Co), vary with the position of the 
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FIG. 8: MAE contributions (upper panel) and orbital magnetic moment variation (lower panel) per atom for Fe4Cos 
and Fe5Co4. Full circles, squares and triangles refer in turn to L z , L x and L y . In the case of Fe-riCos, the Fe is in 
position 1 (denoted by Fei) and Co are in positions: 2, 3 and 4 (denoted in turn by Coi, C02, and C03); in the case 
of FesCo4, the Co is in position 1 and the Fe are in positions 2, 3 and 4. (L Co in Co/Cu(100) is 0.139/^B/atom; L Fe 
in Fe/Cu(100) is 0.08/is/atom). Triangles and full diamonds refer to AEii and AEx, respectively. A-Eii = AE yx and 
AE± — (AE X z + AE yz ) /2. The indices _L and || indicate that the MAE is calculated once for the magnetization switch 
from the normal z towards the in- plane x (y, respectively) direction and then for the in-plane anisotropy between the x 
and the y direction. 



atom within the cluster and increases with increasing coordination number. This example illustrates that the 
preferred magnetization direction of the cluster as a whole results from an interplay of contending tendencies 
of the cluster constituents. 

It is interesting to compare also the configuration averaged MAE, which is thought to be measured in exper- 
imental situations of the 2x2 and 3x3 with those of the pentamer clusters as in Fig. [9l One notices that the 
MAE of the 2x2 and 3x3 clusters have the same tendency as in the monolayer case, whereas for the pentamer 
clusters the situation is different. The larger values for most concentrations in that case are due to the high 
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FIG. 9: Upper panel: MAE per atom for the 2x2, 3x3 and the pentamer cluster, averaged over all configurations with 
respect to the Co concentration. (AE±) = {(AE XZ + AE yz ) /2) corresponds to the average over all configurations at a 
given concentration. The index _L shows that the MAE is calculated with respect to the magnetization along the surface 
normal z and the in-plane x or y direction. Lower panel: anisotropy of the averaged orbital magnetic moment for a 2x2, 
a pentamer and a 3x3 cluster. (AL XZ ) = (L x — L z ), where () refers to the average over all configurations. 



ratio of corner atoms which have a reduced coordination and whose contributions to the MAE are consequently 
the largest. In all the cases we have a change of direction for the easy axis with respect to the concentration. 
For the 2x2 and 3x3 clusters this change occurs at approximately the same concentration (~ 40% Co), while 
for the pentamer cluster, the Co concentration at which the MAE direction changes is higher (~ 70%). In 
the lower panel of Fig. 03 we show that the behavior of the MAE follows closely the averaged orbital moment 
anisotropy. 

The previous sections have focused on the energy difference between the directions along the coordinate axes 
x, y and z only. However, one can expect uniaxial anisotropics of the in-plane magnetization along intermediate 
directions due to the symmetry of a specific cluster configuration. Hence in the following we investigate the 
angular dependence - specified by the azimuth angle ip - of the MAE for selected clusters, in order to illustrate 
the complex behavior of the in-plane anisotropy. 

Fig. \W\ shows the dependence of the total MAE on the azimuth angle for the simple case of one Fe atom 
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FIG. 10: Dependence of the anisotropy energy, AE(ip) = E(ip) — E(ipo), with respect to different directions of the 
magnetization, specified by the azimuth angle ip. The three curves stand for the different configurations (c.f. Fig. I12|l : 
no 1, corner Fe atom (dashed line); no 2, central Fe atom (dash-dotted line), no 3, edge Fe atom (solid line) 



in a 3 x 3 Co cluster. Within this cluster the Fe atom can occupy three inequivalent sites - corner, edge, 
center - and the direction of the easy axis, as well as the MAE strongly depends on the position of the Fe 
atom. For the cluster with the higher symmetry with respect to the x and y axis (Fe atom at center), the 
value of the anisotropy is very close to zero and its angular variation is only due to the cluster symmetry, which 
has a periodicity of 90° (inset of Fig. [TO]) . If the Fe atom is located at the center of an edge the anisotropy 
energy is two orders of magnitude higher. In case the Fe atom sits on an edge parallel to the i-axis, the easy 
magnetization axis is also along the x direction and the hard axis is perpendicular to it along the y direction. 
The periodicity of the angular variation of the MAE in this case is 180°, and the energy needed to switch the 
magnetization between two easy axis is identical to the amplitude. If the Fe atom is positioned at a corner, 
then the easy axis is along ip=45°. Consequently the angular dependence of the MAE can be approximated 
by: AE(ip) = K(C)sm 2 (ip — (po) , where fo is the direction of the easy magnetization axis, depending on 
the configuration symmetry, and K(C) is an anisotropy constant, which is very sensitive to the specific cluster 
configuration C. 

As an additional illustration, Fig. QT] shows the values of the anisotropy constants K{C) and the angles of 
the easy axis for Fe2Coy. This cluster is distinguished from the previously discussed one by the appearance of 
two additional angles at <po «23° and ipo «63° corresponding to symmetry axes of the clusters. The number at 
the left of each symbol refers to the corresponding cluster configurations as given in Fig. [121 

These simple cases illustrate the general features of more complicated ones, namely that the direction of the 
easy axis is determined from the symmetry of the cluster configuration, and that the size of the anisotropy energy 
depends on the specific distribution of atoms. Clearly this means, that by averaging over a large collection of 
clusters, where the exact individual configurations are not known, it is neither possible to predict the anisotropy 
energy, nor the direction of the easy axis of a single nanocluster. 



13 















9 O 




0.4 




Fe 2 Co 7 








(meV/atom) 


0.3 
0.2 










40 


















0.1 






5 
7 


6 O 


10 o 







i 


11 O 
i 


8 


O 

i i 


i 









23 




45 63 
(deg) 


90 



FIG. 11: The anisotropy constants K(C) for eight different configurations of Fe2Co7 cluster. The specified angles (fio 
refer to the direction of the easy axis of magnetization with respect to the cluster configurations. For each case, the 
corresponding configuration number (c.f. Fig. I12p is indicated next to the symbols. 



IV. CONCLUSION 

In this work we have determined the magnetic moments and the magnetic anisotropy of three different types 
of composite FeCo nanoclusters on a Cu(100) substrate by means of ab-initio calculations. The study has 
concentrated on the effect of composition on the cluster properties, neglecting structural relaxation effects. 
While the atomic potentials are approximated as being spherically symmetric, higher moments of the charge 
density have been taken into account. Clusters of three different geometries have been considered - 2x2, cross- 
like pentamer, and 3x3 - and the positions of the Fe and Co atoms have been varied according to all possible 
arrangements for given "concentrations" . 

As in clusters of a pure material, the size of the spin magnetic moments depends on the position in the 
cluster for both, Fe and Co atoms, and consequently on the number of coordinated atoms as has been noticed 
previously by Mavropoulos 36 . In addition, depending on the atoms' position and on the geometry of the cluster, 
the spin-moments show either no, or a high sensitivity to changes in the local environment, i.e., with respect to 
changes of the number or the position of the foreign atomic species. 

Orbital moments are strongly influenced by the position of an atom. As has been found in many previous 
calculation a 22 i 31 ' 34 i 35 ' 37 i 38 ' 39 i 40 i 41 i 42 ' 43 i 44 low coordinated atoms exhibit a significant increase of their orbital 
moments. In the present study it could be shown, that they also depend on the specific cluster configurations 
and on the concentration of the atomic species. Moreover, low coordinated Co atoms, which have the largest 
moments within the surface plane, also exhibit a strong in-plane anisotropy of their orbital moments. Perimeter 
Fe atoms, in contrast, have strongly enhanced orbital moments along the z-direction, perpendicular to the 
substrate surface. 

By averaging over all possible cluster configurations, the total magnetic moments show a variation with respect 
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FIG. 12: Inequivalent configurations of 3x3 clusters containing one or two atoms, respectivley, of a foreign species. Filled 
(empty) circles represent Fe (Co) atoms and vice versa. 



to the concentration of Co atoms, which is very close to that found in monolayers of Fe^Coi-a; as obtained from 
CPA calculations. 

When evaluating the magnetic anisotropy energy it appears that the orientation of the easy axis depends 
on the cluster symmetry, originating from the distribution of the two different magnetic atoms. Then again, 
the size of the anisotropy energy is highly sensitive to the specific arrangements of atoms. Clusters which have 
a preferred in-plane magnetization direction are found to have an unusually large angle dependent anisotropy 
energy. In these cases the direction of the easy axis will not necessarily be along the crystal axes, but is 
determined by the specific distributions of Fe and Co atoms. 

By plotting the MAE, averaged over all configurations, as a function of concentration of one atomic species, 
it is found that the easy magnetization direction of a collection of clusters can be tuned in a similar manner, as 
in thin films of the same materials,^ by varying the concentration of the constituents. 
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